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Abstract. The vast set of near global and continuous atmospheric mea- 
surements made by the SABER instrument since 2002, including daytime 
and nighttime kinetic temperature (T fc ) from 20 to 105 km, is available to 
the scientific community. The temperature is retrieved from SABER mea- 
surements of the atmospheric 15 //m CO 2 limb emission. This emission sep- 
arates from local thermodynamic equilibrium (LTE) conditions in the rar- 
efied mesosphere and thermosphere, making it necessary to consider the CO 2 
vibrational state non-LTE populations in the retrieval algorithm above 70 km. 

Those populations depend on kinetic parameters describing the rate at which 
energy exchange between atmospheric molecules take place, but some of these 
collisional rates are not well known. We consider current uncertainties in the 
rates of quenching of CC^^) by N 2 , O 2 and O, and the 002 (^ 2 ) vibrational- 
vibrational exchange to estimate their impact on SABER T fc for different at- 
mospheric conditions. The T*, is more sensitive to the uncertainty in the lat- 
ter two and their effects depend on altitude. The T*, combined systematic 
error due to non-LTE kinetic parameters does not exceed ±1.5 K below 95 km 
and ±4-5 K at 100 km for most latitudes and seasons (except for polar sum- 
mer) if the T*; profile does not have pronounced vertical structure. The er- 
ror is ±3K at 80 km, ±6K at 84 km and ±18 K at 100 km under the less favourable 
polar summer conditions. For strong temperature inversion layers, the er- 
rors reach ±3 K at 82 km and ±8 K at 90 km. This particularly affects tide 
amplitude estimates, with errors of up to ±3 K. 
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1. Introduction 

The Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) is 
a broadband radiometer which measures limb radiance emitted from the terrestrial atmo- 
sphere in 10 selected spectral bands ranging from 1.27 to 15 /xm. SABER, launched aboard 
the Thermosphere-Ionosphere-Mesosphere Energetics and Dynamics (TIMED) satellite on 
December 7 th of 2001, orbits the earth in a 74° inclined plane. It scans the limb from the 
surface up to ~350km and then back down to the surface in a 110 second sequence with 
a 2 km field of view and with a very high signal to noise ratio (S/N=l at ~ 130 km for 
the 15 /xm channels). SABER daily provides about 1500 atmospheric radiance profiles per 
channel both during day and night covering from 82°N to 53°S (alternating with 53°N to 
82°S every two months). A more detailed description of SABER performance is given by 
Russell et al. [1999]. 

One of the main goals of the SABER experiment is to retrieve kinetic temperature (T fc ) 
from 20 to 105 km with a high accuracy. The temperature is needed not only to study the 
atmospheric dynamics, chemistry and energetics but also to retrieve almost all SABER 
products. In order to derive T*,, two channels are devoted to register CO 2 infrared emission 
at 15 /xm. These are the narrow and the wide band 15 /xm channels (channels 1 and 3), 
with bandpasses from 635 to 710 cm -1 and 570 to 780 cm -1 , respectively. They measure 
over the spectral region of the CO 2 vibration bending mode energy ( 1/2 —667 cm -1 ) with 
noise equivalent radiances (NER) of 2.45 x 10 -4 Wsr -1 m -2 and 3.3 x 10 -4 Wsr -1 m -2 , 
respectively. 
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In a first step, the pressure and the kinetic temperature are retrieved using the two- 

channel technique suggested by Gille and House [1971] assuming LTE (local thermody- 
namic equilibrium). However, since the assumption of LTE breaks down in the meso- 
sphere (see, for example, Lopez- Puertas and Taylor, 2001), the T fc is derived in a second 
step with the algorithm described in Mertens et al. [2001], which accounts for non-local 
thermodynamic equilibrium (non-LTE) effects. Thus, it is important to understand the 
mechanisms that populate the CO 2 vibrational levels emitting at 15 pm. That applies, in 
particular, to the rates at which energy exchanges between atmospheric molecules take 
place. Some of these rate coefficients are still not determined with high accuracy and their 
current uncertainty is the main source of error in the calculation of the level populations. 
That induces non-negligible errors in the retrieved parameters. The magnitude of these 
errors depends on the atmospheric conditions, being larger when the non-LTE effects are 
more important, i.e., under polar summer conditions or for large amplitude temperature 
inversion layers. 

In a companion paper, Remsberg et al. [2008] validate SABER version 1.07 (vl.07 
hereafter) temperature from 20 to 100 km and assess the errors related to the retrieval 
scheme and to the noise and the pointing jitter. They review in detail the errors up 
to the lower mesosphere due to the approximations assumed in the forward radiance 
calculations, uncertainties in the channel 3 radiance calibration, selection of the region 
of pressure registration, and uncertainties in the ozone contribution correction for the 
channel 3 radiance. The overall (root-sum-squared) estimated random and systematic 
temperature errors are respectively 0.3 and 1.4 K at 100 mb, 0.3 and 0.7 K at 10 mb, 0.6 
and 1.9 K at lmb, and 0.7 and 1.5 K at 0.1 mb (see their Table 1). They also provide an 
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estimate of the errors due to the assumed CO 2 profile in the upper mesosphere and lower 
thermosphere (1.3 K at 80km, 3.6 K at 90km, and 1.4 K at 100km), and the noise error 
for a single scan (1.8 and 2.7 K at 80 km, 3.6 K and 8.9 K at 90 km, and 6.7 and 8.9 K at 
100km for mid-latitude and polar conditions respectively). 

The work presented here complements the Remsberg et al. paper in that it assesses the 
errors induced by the current uncertainties in the non-LTE model collisional rates on the 
temperature in the Mesosphere and Lower Thermosphere (MLT), which are added to the 
systematic errors (see their Table 2). We estimate those errors using SABER measure- 
ments taken during typical and extreme conditions, and different seasons. In Section 2 
we describe the main features of the atmospheric emission at 15 //m in the MLT region, 
summarize the SABER non-LTE T* retrieval code and refer to the SABER data used for 
this analysis. In Section 3 we review the collisional rate coefficients used in the non-LTE 
model whose uncertainties most affect the T fc retrieval and report the systematic errors 
on the retrieved T* caused by the current uncertainties in those collisional parameters. 
In the last section, we summarize our results and draw our conclusions. 
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2. SABER non-LTE Kinetic Temperature 

Mesospheric and lower-thermospheric kinetic temperature is retrieved from SABER 
atmospheric 15 /im radiance measurements considering non-LTE effects. Mertens et al 
[2001, 2002] describe in detail the non-LTE retrieval algorithm used operationally to derive 
that kinetic temperature (non-LTE T*, hereafter) above 40 km. This code is composed of 
two main components: (i) the forward model; and (ii) the inversion model. 

The forward calculations at 15 pm are performed with the BANDPAK algorithm [Mar- 
shall et al, 1994] after the populations (i.e., the vibrational temperatures) of the CO 2 
vibrational levels are supplied by a non-LTE model. This model includes all known 
excitation and relaxation mechanisms of the states, i.e., collisional processes (thermal col- 
lisions, vibrational energy transfer, electronic to vibrational energy transfer) and radia- 
tive processes (energy exchange between atmospheric layers, absorption of solar radiation, 
spontaneous and induced emission), as described in Lopez- Puertas and Taylor [2001]. The 
vibrational temperatures (T v hereafter) calculated with the SABER non-LTE model were 
validated against those from the ARC [Winter steiner et al., 1992], the IAA [Lopez- Puertas 
et al, 1986, 1998] and the ALI [Kutepov et al, 2006] non-LTE models from 40 to 160km. 
SABER operational non-LTE and IAA codes solve the radiative transfer problem using 
the Curtis matrix method, whereas the ARC and the ALI (Accelerated Lambda Iteration) 
codes use an iteration technique. All models included the same collisional rate constants 
in the intercomparison. The differences between the models are smaller than 1 K for all 
levels included at all altitudes considered. 

The inversion code consists of two nested relaxation loops. In the inner loop, kinetic 
temperature is retrieved with an onion peeling procedure using a Levenberg-Marquardt 
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approach. This module is iteratively run until the convergence criterion is met, i.e., until 
the difference in the simulated radiance at every tangent height between two successive 
iterations is less than the noise error. In the outer loop the pressure is hydrostatically 
rebuilt and the vibrational temperatures are again calculated. The outer loop convergence 
criterion requires the root-sum-square (RSS) of the difference of two successive T*, profiles 
over all retrieved altitudes be less than the RSS of the expected solution errors (noise- 
induced) over all retrieved altitudes. 

In order to understand the response of the non-LTE retrieval to changes in the collisional 
parameters, we briefly describe here the behavior of the CO 2 vibrational level populations 
and their contributions to SABER 15 /xm channels. We will use hereafter the Herzberg 
notation to refer to the vibrational levels, i.e., (v\,v l 2 ,vz), where vi,V 2 ,vs designate the level 
excited in the symmetric stretching, the bending and the asymmetric stretching modes 
respectively and l the vibrational angular momentum quantum number. We will also 
denote the CO 2 isotopes with the last digit of the atomic weight of each atom, e.g., 626 
isotope for 12 C 16 02 - 

The vibrational levels and bands which, directly or indirectly, affect the CO 2 emission 
at 15 /xm measured by SABER are shown in Figure 1. The fractional contributions at 
non-LTE altitudes to the total SABER 15 /xm channel 1 radiance of the CO 2 vibrational 
bands (see bands included in the model in Mertens et al., 2002) for a standard case are 
shown in Fig. 2. We have estimated those from simulated profiles using the US Standard 
76 atmosphere [Fleming et al., 1990]. Although the fractional contributions depend on 
the atmospheric conditions, these results are representative of the general behavior. The 
626-CO2 V 2 fundamental band (01 1 0-00°0) is the main contributor to the 15 /xm channel 
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total radiance at all altitudes although the first hot band (10°0-01°0, 02 2 0-01 1 0 and 02°0- 

01°0) and the minor isotope (636, 628, 627) fundamental band contributions reach 30% 
at 85 and 75 km, respectively. The contribution from the remaining bands (rest of CO2 
hot bands and ozone v 2 fundamental band) is less than 5% above 70 km. 

SABER temperatures from version 1.07 were retrieved assuming: a) CO2 volume mixing 
ratio (vmr hereafter) from monthly and diurnal averages from the Whole Atmosphere 
Community Climate Model (WACCM) [Garcia et al., 2007 and references therein], b) N 2 , 
O2, nighttime O and daytime O (above 95km) vmr from NRL-MSISE-00 [Picone et al., 
2002]; c) SABER daytime O (below 95 km) as inferred from the O3 retrieval of SABER 
1.27 /am channel [Mlynczak et al., 1993, 2007]; d) the nominal values of the collisional 
rate coefficients (see Lopez- Puertas et al. [2001, Table 6.2]; Kutepov et al. [2006]; and 
discussion below). Also, the non-LTE T fc is merged with the LTE T fc around ~65km. 

Figure 3 shows zonal mean (latitude cross-section) maps of SABER vl.07 kinetic tem- 
perature for equinox (March 18 th 2004) and solstice (July 15 th 2004). We selected these 
two days because yaw maneuvers were performed then, hence, SABER latitudinal coverage 
is the largest, from 82°S to 82°N. Whereas the kinetic temperature in the MLT during 
the equinox is almost latitude-independent (except for the stronger tidal signatures at 
low latitudes and perhaps a slightly warmer lower mesosphere in the northern most lat- 
itudes), the northern summer solstice distribution shows the typical mesopause positive 
temperature gradient towards the winter polar region during the solstice, together with a 
lower mesopause in the summer polar region. Since the atomic oxygen is very important 
in the SABER T*, determination in the upper mesosphere and the lower thermosphere 
(see Sec. 3.2 for a detailed discussion), we also show the zonal mean daytime and night- 
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time atomic oxygen mixing ratio used in vl.07 T*, retrieval for the selected days (Fig. 4). 
Remsberg et al. [2008] shows examples of the CO 2 vmr profiles used. 

Figure 5 shows vl.07 averaged kinetic temperature profiles together with the corre- 
sponding vibrational temperatures for the main CO 2 Vi levels calculated with the SABER 
operational non-LTE model for a typical atmospheric scenario (mid-latitudes) and two ex- 
treme cases (polar summer and polar winter). We calculated those zonal means using all 
SABER temperature profiles measured within 10° latitude bins during January 15 th 2003 
for the mid-latitude conditions and July 15 th 2004 for the polar conditions. The average 
latitude of these zonal mean profiles is 49°S, 79°N and 69°S respectively. 

The behavior of the v<i level vibrational temperatures is mainly driven by the thermal 
structure [Lopez- Puertas and Taylor, 2001]. In the polar summer, the populations of the 
15 pm levels are further away from LTE. In particular, the 626-01 1 0 level (responsible 
for the main contribution to SABER 15 pm radiance) is significantly overpopulated with 
respect to LTE conditions (its vibrational temperature is larger than the T*,) in the polar 
summer from 75 km to 90 km, whereas it is similarly or less populated (its vibrational 
temperature is similar or smaller than the T fc ) at mid-latitudes and in polar winter, where 
the radiative equilibrium conditions are closer to the thermal balance. That happens 
because, besides collisions with other molecules (not very frequent in this low pressure 
region), absorption of radiation emitted at lower warm atmospheric layers increases the 
population of the vibrational levels in that altitude range. In polar summer, the warmer 
stratopause together with the colder mesopause make the vibrational temperatures larger 
than the kinetic temperature around the mesopause. In mid-latitudes and polar winter, 
the mesopause temperature is not as low and, thus, the absorption of upwelling radiation 
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is not enough for the vibrational temperatures to overcome the local kinetic temperature. 

Above about 90 km, the 626-01 1 0 level population is smaller than in LTE under all condi- 
tions because collisions are not frequent enough to counteract the emission to space in the 
01 1 0-00°0 band. Besides, the populations of the 626-01 1 0, the 626-020 (group composed 
by the 10°0, 02 2 0, and 02 2 0 levels) and the 636-01 1 0 levels, depart from LTE at lower 
altitudes (70-75 km) in the polar regions. As we show in Sec. 3, the larger the deviation 
from LTE conditions, the larger the impact of perturbations of the collisional rates and, 
therefore, the larger the error on the retrieved kinetic temperature. 
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3. Kinetic Parameter Uncertainty Ranges and Errors on SABER T*. 

The full set of collisional processes used in SABER CO2 non-LTE model is described 
in Lopez-Puertas et al. [2001]. We have considered here the four collisional non-LTE 
processes that most affect the CC^fy^) populations and hence the retrieved T fc . They are 
summarized in Table 1 together with their nominal rate coefficients (as used in version 
1.07) and the uncertainties assumed in this work. The selection of these values is based on 
considerations made by Winter steiner et al. [1992], Lopez-Puertas et al. [1998], Lopez- 
Puertas and Taylor [2001] and Kutepov et al. [2006] and are discussed below. There 
are other non-LTE processes included in the SABER non-LTE model that have a non- 
negligible uncertainty but they do not yield significant errors in the retrieved temperatures. 

To estimate the T*, errors caused by those uncertainties, we have perturbed the nomi- 
nal rate constants used in the model and performed T fc retrievals. Comparisons with the 
kinetic temperature retrieved using the nominal values in vl.07 (nominal T fc hereafter) 
allow for an estimate of these systematic errors. The response of the retrieval to pertur- 
bations of the collisional rates is to compensate for the change induced in the vibrational 
temperatures by modifying the kinetic temperature to restore the T v s close to their nom- 
inal values so that SABER 15 pm radiance is again matched. As shown below, the T fc 
response to the considered increase or decrease of the rate coefficients is not symmetri- 
cal about the nominal value because the behavior of the non-LTE populations and the 
radiative transfer are non-linear. Besides, the magnitude of these T fc errors depends on 
the thermal structure, more strongly affecting the profiles retrieved for those conditions 
in which the vibrational levels are farther from LTE, e.g., under polar summer conditions 
or when strong inversion layers are present. 
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We have fulfilled those calculations for the two typical days during 2004 in equinox and 

solstice (March 18 th and July 15 th ) whose nominal zonal mean T*, is shown in Fig. 3. We 
have also analyzed the impact of the collisional rate uncertainties for some extreme cases 
in which we expect to find larger errors, e.g, those with strong inversion layers in the MLT 
or a high stratopause. We have studied the effect of the collisional rates of processes 3 
and 4 independently whereas we have considered the combined effect of processes 1 and 
2, since both produce a similar effect. 

In the next sections, we justify the selection of the uncertainty ranges for each individual 
process shown in Table 1 and analyze in detail their effect on the SABER T*,. Although 
not a non-LTE parameter, the uncertainty in the atomic oxygen also induces T*, errors 
through the non-LTE population calculations. Those errors are also briefly discussed 
below. The errors due to uncertainties in the CO 2 vmr, however, affect the forward 
radiance through the CO 2 density but not through the non-LTE populations since the 
effect of CO 2 abundance on the vibrational temperatures is small [Lopez-Puertas and 
Taylor, 2001]. Those errors are assessed in [Remsberg et al., 2008]. 

3.1. Processes 1 and 2: C0|(v2)+M(N2,02)^C0|(v2— 1)+M(N 2 ,02) 

Thermal collisions of carbon dioxide with molecular nitrogen and molecular oxygen 
(processes 1 and 2) keep the states in LTE up to the middle mesosphere. The values 
for the rate coefficient reported in the literature for quenching of CXfy^) by N 2 are 
summarized in Fig. 6. Taine et al. [1978] and Allen et al. [1980] measured the rate of 
collisions with molecular nitrogen at atmospheric temperatures but only for low energy 
vibrational levels. There are also laboratory measurements of this coefficient reported in 
Simpson et al. [1977], Allen et al. [1977], Taine and Lepoutre [1979], Lunt et al. [1985] 
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and Siddles et al. [1994], All measurements gave similar values within their experimental 
errors. Simpson et al. [1977], Taine et al. [1978] and Siddles et al. [1994] also provided the 
rate of CCU^) relaxation in collisions with molecular oxygen (see Fig. 7). The nominal 
values for both coefficients (&co 2 -air) used here are the parameterizations of Wintersteiner 
et al. [1992] obtained by fitting the aforementioned experimental data. Based on those 
experiments, we have assumed a [— 25%,25%] uncertainty in the rates of both processes. 
Uncertainties in the temperature dependence are also included in that range. 

Figure 8 shows zonal mean maps of the errors in the retrieved kinetic temperature due to 
this process for the two typical days during equinox (left panels) and solstice (right panels) 
studied here. As with the other collisional rates, the operational T*, retrieval responds as 
expected to the perturbations of the rate constants: the maps for equinox and solstice at 
the top of the figure are almost completely anticorrelated with the corresponding maps 
at the bottom. 

The V 2 fundamental band upper vibrational levels (01 1 0) of all the CO 2 isotopes and the 
first hot band upper levels, main emitters at 15 pm in the middle and upper mesosphere, 
remain close to LTE up to ~85-90km under typical mid-latitude conditions (see Fig. 5-a). 
The uncertainty range of the rate of 002 (^ 2 ) quenching by N 2 and O 2 is not large enough 
to change the altitude of the LTE breakdown significantly and, therefore, the effect on the 
SABER retrieved T*, is small. That is also the case for polar winter and at all latitudes 
during the equinox (see Fig. 8), where the errors are always smaller than ±0.5 K. 

In polar summer, the 01 1 0 levels start to depart from LTE at lower altitudes, in particu- 
lar those of the minor isotopes. Hence, the retrieved kinetic temperature is more affected 
than in mid-latitudes. The uncertainty of these rates produces maximum errors of ±2 K 


DRAFT 


September 18, 2008, 9:40am 


DRAFT 



X - 14 GARCIA-COMAS ET AL.: SABER TK ERROR DUE TO NON-LTE PARAMETERS 

in polar summer conditions at 85 km. It is worth pointing out that, even though the errors 

are small at most latitudes, the effect of increasing or decreasing the collisional rate is 
the opposite for the summer and the winter scenarios below 90 km. This is because the 
626-01 1 0 state is more populated than in LTE in the polar summer mesosphere while less 
populated in polar winter (see Figs. 5-b and c). 

3.2. Process 3: C0*(v 2 )+0( 3 P)^C0*(v 2 -l)+0( 3 P) 

Process 3 plays an important role in the radiative cooling of the terrestrial thermosphere 
since it exerts a major influence on the CO 2 CU 2 ) state populations in this region. Until the 
beginning of the past decade and due to the lack of direct measurements, there were no re- 
liable data for this collisional coefficient. However, several values of this rate have been re- 
ported in the literature since then, some deduced from atmospheric observations (pointing 
to the largest values) and others measured in the laboratory (suggesting the lowest values) 
(see Fig. 9). Values at 300 K range from 6xl0 -12 cm 3 s -1 [Sharma and Wintersteiner 1990] 
(deduced from the SPIRE rocket atmospheric observations) to 1.5xl0 -12 cm 3 s -1 [Shved 
et al., 1991, and Pollock et al., 1993] (from laboratory measurements). 

Khvorostovskaya et al. [2002] carried out the first laboratory measurements in the 206- 
358 K range reporting 1.5xl0 -12 cm 3 s -1 . More recently, Castle et al. [2006] have obtained 
a slightly larger value (1.8x 10 -12 cm 3 s -1 ) from measurements taken at room temperatures. 
As an exception from other values found in the literature, Lilenfeld [1994] measured a very 
low rate of 0.5xl0 -12 cm 3 s -1 . Lopez-Puertas et al. [1992] deduced a value at 300 K of 3- 
6xl0 _12 cm 3 s _1 from ATMOS measurements and Grossmann et al. [2002] corroborated 
the upper value of Sharma and Wintersteiner using CRISTA-2 measurements (after com- 
parison with MSIS data), i.e., both closer to the upper limit previously reported. Gusev et 
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al. [2006] also claimed that the Sharma and Wintersteiner rate constant used to retrieve 
kinetic temperature from CRISTA-1 and CRISTA-2 measurements was consistent with 
the global CRISTA 15 fim measurements above 120 km. The lowest value derived from 
atmospheric measurements so far reported is that by Vollmann and Grossmann [1997], 
who, after the analysis of SISSI data, suggested a rate of 1.5xl0 -12 cm 3 s -1 . However, 
Gusev et al. [2006] pointed out that this low value was most probably due to an unreal- 
istically large CO 2 number density in the lower thermosphere (more than 4 times larger 
than usual) used in the analysis. 

The nominal value used in vl.07 for this rate constant (£002-0) is that reported by 
Sharma and Wintersteiner, corroborated by all other atmospheric measurements, includ- 
ing their suggested temperature dependence. Although the laboratory measurements of 
&C02-0 grouped around a lower rate, the atmospheric measurements do not support 
that value, clearly pointing to a rate of collisions about four times larger, i.e., that of 
Sharma and Wintersteiner. Since vl.07 SABER T*, retrievals use that upper limit rate, 
we have assumed the uncertainty range provided by Sharma and Wintersteiner, i.e., [- 
50%, +50%]. The uncertainty of the temperature dependence of this rate is also included. 

The quenching of CO 2 by atomic oxygen is a key process for determining the kinetic 
temperature from measurements at 15 /am in the upper mesosphere and lower thermo- 
sphere, where O is abundant. The latitudinal and seasonal dependence of the errors due 
to the kco 2 -o uncertainty is shown in Fig. 10. In order to explain these errors, it is useful 
to understand first the response of the vibrational temperatures to kco 2 -o perturbations. 

Lowering the rate of collisions between C02(f2) and atomic oxygen leads to level 
populations farther away from LTE. Figure 11-a shows the effect on the vibrational tem- 
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peratures of the main CC^^) vibrational states when dividing the rate by two for a 

typical mid-latitude scenario. The altitude at which the C02(t>2=l) levels of all isotopes 
depart from LTE barely changes because the effect on their populations is only important 
a few kilometers above the altitude of LTE breakdown. Since the T v s of the 626-01 1 0 and 
636-01 1 0 states are smaller than the T*,, reducing the coupling to the kinetic temperature 
by lowering the rate kco 2 -o further decreases them above 80-85 km (by 3 K and 1 K re- 
spectively at 95 km). Although the opposite happens to the 01 1 0 levels of the 628 isotope 
(its T v increases 2-3 K at 95 km), the overall effect is to decrease the 15 /im forward limb 
radiance. This is compensated by an increase in the retrieved T*, to match the SABER 
radiance, i.e., the induced error is positive. The errors in the kinetic temperatures in the 
positive sense are 1-2 K from 84 to 96 km, 5K at 100 km and up to 11 K at 104 km for 
most latitudes except for polar summer. 

In contrast, increasing kco 2 -o brings the 01 1 0 level populations closer to LTE. The T v s 
of all levels increase, except for the weaker 628 and 627 isotopic 01*0 levels that decrease 
between 80 and 100 km. In that case, the resulting T*, decreases at all altitudes with 
respect to the nominal case. The errors in the negative sense are smaller (in absolute 
value) than IK below 96km, and from -1 to -2.5 K from 100 to 104km (see Table 2). 

Except for polar summer conditions (discussed below), the response of the retrieved T& 
to changes in this rate is similar for all latitudes and seasons (see Fig. 10). For polar winter 
scenarios, the effect appears a few kilometers lower since the 01 1 0 states depart from LTE 
at lower altitudes but, as in mid-latitudes, increasing kco 2 -o leads to a decrease in the 
retrieved T fc at all altitudes (see also Table 4). Also, as we discuss in more detail later, 
the errors are larger in the tropics mainly due to tidally-induced variability in the kinetic 
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temperature vertical structure. That is also the case for the three maxima in the T*, error 
(3-4 K) found for equinox conditions at 85-90 km around 70°S, 10°S and 50°N (see Fig. 10, 
top- left), not typical of this season but caused by inversion layers at those latitudes (see 
Pig. 3). 

For polar summer conditions, the levels are overpopulated with respect to LTE 
around the mesopause region and less populated at altitudes above around 95-100 km 
(see Fig. 11-b). Therefore, when decreasing the kco 2 -o rate by 50%, the vibrational 
temperatures of all 15 fim CO 2 levels increase in the altitude range of 80 to ~95 km (from 
IK for the 626-01 1 0 level at 85km to 6K for the 01 1 0 minor isotopic levels at 92km) 
but decrease in the region above (up to 12 K at 100 km for 626-01 1 0). Accordingly, the 
retrieved kinetic temperature decreases with respect to the nominal T*, in the mesopause 
region (-2.6 K at 84 km) and increases above (from 6.6 K at 92 km to 40 K at 104 km) (see 
also Table 3). The opposite happens when increasing the collisional rate by 50%, which 
yields errors of IK at 84km and from -1.6 K at 92km to — 10 K at 104km respectively. 
The retrieved T*, is almost insensitive to this rate coefficient around 90 km and also below 
80 km (the latter because atomic oxygen abundance rapidly decreases). 

The average global SABER T fc errors due to process 2 uncertainty are similar to those 
reported for other experiments, e.g., CRISTA [Gusev et al., 2006], being slightly smaller 
below 90 km and larger above. 

In this work, we have assumed a kco 2 -o of 6x 10 -12 cm 3 s -1 at 300 K with an uncertainty 
of [-50%, +50%]. All determinations of this rate constant from atmospheric measurements 
lie within that range. However, laboratory measurements point to a value around half of 
the lower limit considered here. Even if the aeronomic evidence does not support that 
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low value (see above), we have also calculated the change in the kinetic temperature that 

would result when lowering the nominal rate constant by a factor of 4. The differences 
in T k are then two to three times larger than those obtained when dividing &C02-0 by 2 
(see Tables 2-4). 

Although, the atomic oxygen number density is not a non-LTE parameter, the sensitiv- 
ity of T*; to its uncertainty is similar to that of the uncertainty in the CO2-O quenching 
rate. As mentioned in the previous section, vl.07 T fc retrieval assumes daytime 0 vmr 
below about 95 km derived using SABER O3 retrieval of 1.27 /xm measurements and the 
0 2 ( 1 A) model of Mlynczak et al. [1993, 2007] and NRL-MSISE-00 atomic oxygen vmr 
everywhere else. MSIS atomic oxygen vmr is modeled from measurements of the total 
oxygen amount ([O] + 2[0 2 ]) in the atmosphere ( Picone et al., 2002). Picone et al .' s 
comparison of NRL-MSISE-00 molecular oxygen with that of MSIS-90 shows differences 
of ~5%, which translate to a ~10% difference in atomic oxygen between the two versions. 
Gumbel [1997] comparisons of MSIS-90 atomic oxygen with rocket-borne measurements 
for several seasons and latitudes show a ~40% difference below 95 km for nighttime and 
±30-40% above. More recently, Mlynczak et al. (2004) radiance simulations at 63 /xm us- 
ing MSIS atomic oxygen agreed within 10% with FIRS-2 measurements. Taking all those 
considerations together, we can assume an uncertainty in the atomic oxygen abundance 
used in SABER T& retrieval of ±50%. The magnitude of those errors is thus similar to 
that produced by the kco 2 -o rate. This error is additional to those due to the kinetic 
parameters. 
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3.3. Process 4: CO*(u 2 )+CO^(u 2 ) ^ CO|(v 2 — l)+CO i(v' 2 + 1) 

The set of v 2 vibrational energy exchanges between C0 2 molecules (process 4) is impor- 
tant for the redistribution of the u 2 quanta among C0 2 molecules. The most important 
V-V process considered in this set is the exchange among C0 2 isotopes (i> 2 =l and v' 2 —0, 
and i^j in process 4 of Tablet), in particular, when one of the colliders is the main iso- 
tope. It also includes the most likely V-V exchange involving the higher excited levels 
V 2—2 and p 2 =3, i.e., that with the lowest energy level of the same isotope (i4=0 and i—j ). 
In this second sub-set of processes, the most important are the V-V exchanges within 
the main isotope. Process 4 mainly affects the populations of the 01 1 0 levels of the less 
abundant isotopes and, to a lesser extent, that of the 020-group (composed of the coupled 
triad 10°0, 02 2 0, 02°0; see Fig. 1) and higher levels of the main isotope from 75 to 90 km. 
Above that altitude, quenching by atomic oxygen dominates. 

In an early experiment, Rhodes et al. [1968] indicated a rate of relaxation of the 10°0 
level in V-V C0 2 exchange of 1.15xl0 -11 cm 3 s -1 , although they pointed out that they 
could not distinguish if the 02 2 0 level was involved. Later interpretations of that exper- 
iment [ Seeber , 1971] proposed that the rate measured by Rhodes et al. was more likely 
not related to our process 4. Additionally, Seeber calculated a rate for the 02 2 0 state 
relaxation leading to a combined relaxation rate of the 020 states of 1.6xl0 -11 cm 3 s -1 . 

More recently, the rate at which the 020 triad levels relax to the 01 1 0 level through V-V 
exchange was measured by Huddleston and Weitz [1981] , who reported an overall effective 
rate of 1.2xl0 -11 cm 3 s -1 . Dang et al. [1983] measured, instead, the rate for the relaxation 
of the individual 020 levels at 280 K, concluding that the relaxation of the 02 2 0 level 
(4.7xl0 _11 cm 3 s _1 ) is about 20 times faster than those of the other two levels of the triad 
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(10°0 and 02°0). Assuming the levels in the triad are in LTE among themselves (which 

is a good assumption since their collisions with N 2 and O 2 which couple these levels are 
about two orders of magnitude more frequent than the V-V exchange), the rates of Dang 
et al. yield an effective rate for the relaxation of the 020 triad of 2.4xl0 -11 cm 3 s -1 , i.e., 
twice the value measured by Huddleston and Weitz. Based on theoretical arguments, Orr 
and Smith [1987] agreed with the laboratory measurements of Dang et al. 

Considering the harmonic oscillator law, the rate for the V-V exchange of the CO 2 
isotopes in the fundamental V 2 —I state can be approximated by half of the 02 2 0-01 1 0 
relaxation rate (i.e., 2.35xl0 -11 cm 3 s -1 according to the measurements of Dang et al.). 
Therefore, the rate at which the first excited level 01 1 0 relaxes in V-V collisions is almost 
equal to that at which a molecule in the 020 triad relaxes. 

Kutepov et al., [2006] showed that kinetic temperature measurements of falling spheres 
are consistent with SABER vl.07 T fc , retrieved setting the V 2 exchange rate among the 
CO 2 isotopes (process 4 with V 2 —I and v' 2 —0) to 2.4xl0 -11 cm 3 s -1 and the rate constant 
for the effective relaxation of the 020 triad to 1.2xl0 -11 cm 3 s -1 , i.e., two times smaller 
than for the inter-isotope exchange. Those are the nominal rate constants used here. This 
combination represents an intermediate case between that provided by Dang et al. and 
that by Huddleston and Weitz, actually slightly closer to that of Dang et al. since the 
inter-isotope V 2 —I exchange has a greater effect. We have estimated the errors assuming 
the same rate constant for both the isotope relaxation and the 020 triad relaxation in V-V 
i >2 collisions (denoted by few), setting it to 1.2xl0 -11 cm 3 s -1 (based on Huddleston and 
Weitz) and 2.4xl0 -11 cm 3 s -1 (based on Dang et al). 
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The CO 2 V-V V 2 exchange mainly affects the population of the minor isotope V 2 states 
and, to a lesser extent, the 020 and higher V 2 levels of the main isotope in the upper 
mesosphere. Those levels are significantly more affected than the 01 1 0 state of the main 
isotope because the latter contains the majority (94%) of the V 2 excited molecules and, 
thus, this process produces only a small relative change of its population. The inter- 
isotope V-V exchange act to depopulate the minor isotope states in favor of the 01 1 0 
state of the main isotope. In case of V-V exchange within the same isotope, this process 
depopulates the higher energy level involved. 

The effect is largest under polar summer conditions because the levels are then farther 
away from LTE. Therefore, as pointed out by Kutepov et al. [2006], it is important to 
include this process in the retrieval to obtain an accurate kinetic temperature in the polar 
summer mesopause. In Fig. 12, we show the effect of changing few on the vibrational tem- 
peratures for a typical polar summer temperature profile. When considering the smallest 
fc vv rate constant for all V 2 exchange processes(i.e., reducing only the inter-isotope V-V 
exchange rate by a factor of two but keeping the nominal rate for V-V exchange within 
the same isotope; see discussion above), only the vibrational temperatures of the 01 1 0 
levels of the weakest isotopes are significantly affected. The populations of these levels, 
mainly driven by absorption of radiation from the warmer region below, tend then to in- 
crease at altitudes from 75 to 90 km (where their contribution to SABER 15 pm radiance 
is important) since they become less coupled to the main isotope 01 1 0 state, which has a 
lower vibrational temperature. The T v decrease in the example shown is 4 K for 636-0 1*0 
and up to 8 K for 628-01 1 0. Therefore, the retrieved T*, for a slower few rate is smaller 
than the nominal temperature at those altitudes. 
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The errors in the kinetic temperature for the equinox and the solstice days considered 

here (March 18 th and July 15 th , 2004) are shown in Fig. 13 (top panels). They do not 
exceed 0.5 K at any altitude, latitude or season except close to the polar summer, where 
the T*; error around 84 km ranges from -2K at 40° N to -7K at 82° N (during July 15 th 
2004). Although this process barely affects the vibrational level population above 90 km 
(see Fig. 12), the temperature slightly increases above 90 km around the polar summer 
(by 2 K) for a slower rate. The temperature decrease around the mesopause resulting 
from the decrease in the kyy rate constant causes a decrease in pressure which, through 
hydrostatics, extends to higher altitudes. This decrease in pressure above 90 km is then 
compensated by a small increase in temperature in order to restore the total density (and 
thus CO 2 density) to match the measured radiance. 

When considering the largest k vv rate constant (2.4x10 -11 cto 3 s -1 ) for all exchange 
processes (i.e., doubling the fc vv involving the v%>2 states and keeping the nominal fc vv for 
the inter-isotope exchange; see discussion above), only the population of the 020-group 
and upper levels are significantly affected (see Fig. 12). The vibrational temperature of 
the 020 levels tends to decrease (2K at 85km in polar summer conditions). The effect 
on the retrieved T fc is negligible under most atmospheric conditions except at latitudes 
close to the summer pole (see Fig. 13, bottom), where T*, increases at 80-90 km causing 
the largest error of 3 K at 84 km. The effect of doubling the rate of the process involving 
the 020 and higher states is smaller than halving the inter-isotope V-V V 2 —I exchange 
rate because the latter is a more important process. 
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3.4. Overall errors 

Tables 2, 3, and 4 summarize the errors on SABER T*, due to the uncertainty (as 
indicated in each column) of each non-LTE kinetic parameter for mid-latitudes, polar 
summer and polar winter, respectively. The mid-latitude errors are an average of the errors 
estimated for T*, profiles measured at latitudes within [35°N,55°N] and [35°S,55°S] during 
equinox and solstice (March 18 th and July 15 th 2004). The errors for polar conditions are 
the average of the errors of profiles within 10° around the corresponding polar latitude in 
the solstice case. For each rate constant, the column labeled with ‘typ’ shows the typical 
error, estimated with the quadratic average of the errors in both directions. 

We also report in those tables the extreme errors (max - and max + ) obtained by setting 
the rate constants of the four processes 1-4 within their uncertainty range to produce the 
largest deviations (in both senses) from the nominal case. Those errors represent the 
less favourable and rather conservative cases because all rate constants used then act on 
the nominal T*, in the same direction. The set of constants used for each case to obtain 
max - and max + are also indicated in each table. In mid-latitudes and polar winter, 
the maximum deviation from the nominal temperature profile max - is produced for the 
maximum fcco 2 -air and &C 02-0 values (within the uncertainty ranges provided here), and 
the minimum &w value below 88 km and its maximum value above that altitude. The 
maximum deviation in the opposite sense (max + ) is obtained by setting &co2-air an d 
&co 2 — o to their minimum values, and fc vv to its maximum value below 88 km and its 
minimum value above. For polar summer, the maximum deviation from the nominal 
kinetic temperature (max - ) is obtained with minimum fc<x>2-air, &C02-0 and k vv values 
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below 88 km and maximum values above, and vice versa for the maximum deviation in 

the opposite sense (max + ). 

The combined error of all rates, i.e., the root-sum-square of the typical errors due to 
each parameter uncertainty, which represent the (typical size) estimated total errors in 
SABER vl.07 kinetic temperature due to uncertainties in the non-LTE parameters is also 
listed in Tables 2, 3, and 4. The largest typical total errors correspond to polar summer, 
where the most significant occur around the mesopause (±6 K) and above 92 km (>6 K, in 
absolute value). At mid-latitudes, the typical errors are smaller (in absolute value) than 
±1.5 K below 96km and smaller than ±8K below 104km. In polar winter, the errors are 
slightly larger than at mid-latitudes above 96 km, where they reach ±5 K and ±7 K at 
100 km and 104 km respectively. 

It is worth noting that the change of the T fc errors with the solar zenith angle (SZA) 
is small. On the one hand, even if the CO 2 upper levels are indirectly pumped by solar 
radiation, the 626-01 *0 level (main responsible for SABER 15 /xm radiance; see Fig.2) is 
barely affected by the SZA [Lopez-Puertas and Taylor, 2001]. On the other hand, the 
variation of atomic oxygen with SZA is largest below 85-90 km, where the effect of the 
rate uncertainties is small (except for polar summer around 85 km). As a consequence, the 
dependence of the errors on SZA is a second order effect as compared with the dependence 
on the temperature structure, i.e., latitude or strong inversion layers. The differences in 
the overall error between maximum and minimum SZAs for a similar temperature profile 
are smaller than 0.5 K and 1 K for mid-latitudes and polar summer, respectively. 
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3.5. Errors in Tide and Inversion Layer Temperature Amplitudes and in other 
Particular Cases 

In the previous sections we have addressed the effect of the kinetic parameters on the 
retrieved kinetic temperature for scenarios with smooth temperature profiles. However, 
these effects strongly depend on the actual vertical temperature structure and, hence, 
on the presence or absence of inversion layers. This applies in particular to the regions 
where the tides are stronger, i.e., at low latitudes, or when some other heating or cooling 
mechanisms produce oscillating signatures in the temperature profile. 

A typical example for these conditions is that of Fig. 14, where we show the SABER 
kinetic temperature together with the vibrational temperatures for the main CO 2 v% levels 
for an individual profile close to the equator. The kinetic temperature abruptly increases 
around 85 km by 70 K in only 5 km. Under these conditions, the quenching of 002 (^ 2 ) by 
other molecules is not fast enough to bring their populations to LTE, resulting in a large 
departure from LTE and, hence, the vibrational level populations are more sensitive to the 
uncertainties in the non- LTE collisional rates. The maximum errors in T*, are found at the 
altitude of crests and troughs of the inversion layer. Their magnitude depends not only 
on the inversion layer amplitude but also on its altitude. For this individual example, the 
maximum errors in the kinetic temperature caused by the uncertainty in the quenching 
rate of 002 (^ 2 ) by N 2 and O 2 (&co 2 -air)> by O (&co 2 -o)> and the V-V ^-quanta exchange 
rate (fc vv ) are, respectively, ±1.4 K (at 82 km), ±8K (at 90km) and ±2.8 K (at 82 km) 
(see Fig. 15). The overall errors in the crests and troughs of this particular profile are 
±3.5 K at 82 km, ±8K at 90 km and ±1.6 at 100 km. 

These enhanced T fc errors around the inversion layers might affect the tidal analysis 
based on SABER temperatures. We have estimated the errors affecting diurnal tide 
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amplitudes from the differences between the mean temperature of the ascending segment 

scans of an orbit minus that of the descending one at equatorial latitudes. We performed 
these calculations for SABER measurements taken during the north viewing phase of 
July 15 th 2004. The average temperatures for the ascending and descending portions are 
shown in the leftmost panel of Fig. 16. For this particular period, the equatorial local 
time for the ascending scans was 11 p.m. whereas that for the descending ones was 2 p.m., 
i.e., about 9 hours apart. Figure 16 also shows the ascending minus the descending mean 
temperatures divided by two (referred hereafter as amplitude) together with the overall 
effect of the fcco 2 -air, &C 02-0 and Kv rate uncertainties (center panel). The oscillating 
structure of the difference, with a vertical wavelength of about 20 km, and with increasing 
amplitude in the mesosphere is typical of a diurnal tide signature. The error in the 
amplitude is shown in the rightmost panel. The errors in the amplitude are negligible 
below 80 km. For the amplitude of ~16K at 88 km of our example, the overall error is 
±1.5 K, mostly due to the uncertainty in the V-V exchange rate and to a much lesser 
extent to the rate of CO 2 quenching by N 2 and O 2 . The amplitude at 95 km has an error 
of ±3K, which is almost totally due to the uncertainty in the CO 2 -O quenching. The 
errors in the phase are negligible. 

An additional source of error in the estimation of the tidal amplitude might come from 
an inaccurate inclusion of the tide effect on the atomic oxygen used in the T*, retrieval. 
In principle, the effect through the total density (with which the O density is calculated 
from the O vmr profiles) is taken into account because SABER measures pressure and 
temperature. Additionally, the O vmr profiles also include, to a first order, the tidal effects 
since, for daytime cases below 95 km, it is simultaneously measured by SABER and, for 
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any other case, it comes from NRL-MSIS-00, which also incorporates tides [ Picone et 
al, 2002], Nevertheless, we have checked how well the tide shows up in the O vmr for 
our example in Fig. 16. If chemical effects on O are negligible (valid assumption above 
85-90 km [Smith et al., 2003]), the O vmr deviation from the daily average produced 
by a tide can be estimated from the T*, deviation, once the O vertical gradient and 
the static stability are known (see Expression 5 in Smith et al. [2003]). Accordingly, 
the tide of our example, which causes variations of ±16 K at 88 km and 95 km, should 
produce a O vmr variation (in phase with the T fc variation) of ±50-70% and ±30-40%, 
respectively. The deviations of both the ascending and the descending average O vmr 
profiles used in the retrieval with respect to the daily average are 110% at 88 km and 
40% and 95 km, also in phase with the T*, variation. Hence, the effect of the tide on the 
atomic oxygen is correctly represented at 95 km but it is overestimated at 88 km by 50% 
(although within the estimated O uncertainty; see Sec. 3.2). That means that the O vmr 
at 88 km is underestimated by 50% for the descending average and overestimated by 50% 
for the ascending average, which translates into an underestimation of T*, of 1 K and 4 K 
respectively, i.e., both act in the same sense. Consequently, the impact on the amplitude 
of inaccuracies in the inclusion of O vmr tidal effects is an underestimation of 1.5 K at 
88 km but no error at 95 km 

There are other particular scenarios where the thermal structure differs strongly from 
the typical cases studied in previous sections. That is the case for the anomalous high 
stratopause measured by SABER at high northern latitudes during January and February 
of 2004 and 2006. This feature was attributed by Siskind et al. [2007] to a highly disturbed 
lowermost stratosphere which blocked the gravity waves impeding the warming around 
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the usual stratopause altitudes. The errors due to non-LTE rate constants under these 

conditions behave similarly to those around the inversion layers. However, these errors can 
extend downwards depending on the stratopause altitude (as they did with the inversion 
layer altitude). Figure 17 shows the zonal mean kinetic temperature at 80°N for January 
15 th 2004, together with the vibrational temperatures for the main CO 2 v? levels. The 
altitude of the 260 K-stratopause was located at 80 km. The CO 2 v? state populations 
depart strongly from LTE and they do so from uncommonly low altitudes (from 60 km for 
the nominal case). The effect of the uncertainty in the few rate on the kinetic temperature 
(not shown) starts being important at altitudes as low as 60 km, and increases from ±1.3 K 
at 64 km to its maximum of ±2 K at 75 km in this case. It is negligible above 80-85 km. 
The effect of the CO 2 -O quenching rate uncertainty is noticeable above 75 km, reaching 
its maximum of ±5K at 84 km. The effect of the fcco 2 -air rate uncertainty is always 
smaller than ±1K. Thus, the errors under these circumstances become important at 
lower altitudes and are larger than for typical T* profiles. 

In the extreme cases presented in this section, the error in the SABER T* might further 
affect other atmospheric quantity retrievals that need accurate knowledge of the tempera- 
ture profile. That is the case of mesospheric ozone and water vapor retrieved from SABER 
9.6 and 6.3 /im radiance, respectively. The impact of the T* errors in the MLT region on 
other SABER products is outside the scope of this paper and will be subject of a future 
analysis. However, it is important to note that, even considering the enhanced errors 
found under these extreme conditions, the non-LTE SABER T*. retrievals in the MLT are 
more accurate precisely due to the use of non-LTE algorithms. 
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4. Summary and Conclusions 

At the time of this publication, the SABER instrument has provided global middle at- 
mosphere kinetic temperature derived from 15 n m limb radiance measurements for more 
than six years. Non-LTE algorithms are used to retrieve the temperature in the meso- 
sphere and lower thermosphere from SABER measurements [ Mertens et al., 2001]. The 
SABER vl.07 kinetic temperature is validated in Remsberg et al. [2008], who, together 
with the error estimates due to non-LTE parameters presented in this work, provide an 
estimate of the systematic error due to CO 2 abundance uncertainty (1.3 K at 80 km, 3.6 K 
at 90 km and 1.4 K at 100 km) and of the noise error (1.8 K at 80 km, 3.6 K at 90 km and 
6.7 K at 100 km for mid-latitudes, and 2.7 K at 80 km, 8.9 K at 90 km and 8.9 K at 100 km 
for polar summer). In this paper, we have evaluated the systematic errors in the SABER 
operationally retrieved kinetic temperature (version 1.07) due to the current uncertainties 
in the collisional rates of the non-LTE model for the typical mid-latitude conditions and 
also for extreme conditions. The summary of our findings follows: 

1. The response of the operational SABER T fc retrieval scheme to perturbations of the 
collisional rates behaves as expected from non-LTE theory. This is confirmed by the fact 
that increasing the rates produces the opposite response to that when decreasing the rates 
at all altitudes, latitudes and seasons. The systematic error estimations performed in this 
work for typical and extreme atmospheric conditions can, in consequence, be interpolated 
to other intermediate atmospheric conditions. 

2. The overall error in the retrieved kinetic temperature is around ±1-2 K below 95 km 
for mid-latitudes and polar winter, and below 80 km for polar summer. The error is around 
±4K at 100 km in mid-latitudes and ±5K in polar winter. In polar summer, the error 
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is significantly larger, ±6K at 85 km, but decreases to ±1K at 88 km and increases to 

±18 K at 100 km. 

3. Under most conditions, the overall T fc error due to the non-LTE kinetic parameters 
is mainly due to the uncertainty in the quenching rate of 002 (^ 2 ) by atomic oxygen. 
Under polar summer conditions, the uncertainty in the rate of V-V ^-quanta exchange 
also contributes significantly to the overall temperature error below 90 km. 

4. Below 74 km the retrieved T* errors are smaller than 0.5 K in both the typical and 
the extreme atmospheric scenarios considered. 

5. Temperatures retrieved in polar summer scenarios have the largest errors because 
the CO2 vibrational levels populations are farther away from LTE conditions and therefore 
they are more sensitive to the uncertainties in the collisional rates. 

6. Except for situations in which a vertical temperature structure in the mesosphere and 
lower thermosphere is significant, the errors in equinox conditions are latitude independent 
and similar to those at mid-latitudes. 

7. The errors are somewhat larger when the temperature profiles have pronounced 
vertical structure, even outside the polar summer. The overall errors in those cases can 
rise up to ±3 K around 80 km and ±8 K at 90 km if crests or troughs are present at those 
altitudes. This affects the estimate of wave temperature amplitudes from SABER data, 
which can be affected by errors due to non-LTE kinetic parameter uncertainties of ±3 K. 

In order to reduce these systematic errors and raise the altitude at which the SABER T*, 
is accurately derived, it is essential to know with a better accuracy the CO 2 -O quenching 
rate and, in particular, to understand why the measurements in the laboratory differ by 
a factor of 4 with respect to the ones derived from atmospheric observations. It is also 
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necessary to determine the CO 2 V 2 vibrational energy transfer rate with a higher accuracy 
since the published measurements, which date from the eighties, differ by a factor of 2. 

In spite of the error analysis presented here, the kinetic temperature measured by 
SABER provides an unprecedented opportunity to achieve a deep understanding of the 
physics of the mesosphere and lower thermosphere as, for the first time, it is operationally 
derived in that atmospheric region from continuous space flight instrument measurements. 
Additionally, the SABER experiment is providing temperature retrievals that are more 
nearly correct for all of the extreme situations presented in this work thanks to its non-LTE 
algorithm. 
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Figure 1 . Energy of the main CO 2 vibrational states. The arrows show the main CO 2 bands 
emitting at 15, 10 and 4.3 /im (FB: Fundamental band; FH: First hot band; SH: Second hot 
bands; TH: Third hot bands). 

Figure 2. Fractional contribution of CO 2 vibrational bands to SABER Channel 1 (15 /im) 
simulated radiance for a typical mid-latitude profile. (Solid: 626 fundamental band; Dashed: 626 
first hot bands; Dash-dotted: 626 second hot bands; Dotted: 636+628+627 fundamental bands; 
Dash-dot-dotted: Remainder (see text for more details). 

Figure 3. SABER vl.07 kinetic temperatures for equinox (left: March 18 th 2004) and solstice 
(right: July 15 th 2004). The contours are every 10 K. Numbers at bottom x-axis indicate the 
number of averaged profiles. 

Figure 4. Contour plots of daytime (top) and nighttime (bottom) atomic oxygen mixing ratio 
(in ppmv) used in SABER vl.07 kinetic temperature retrieval for equinox (left: March 18 th 2004) 
and solstice (right: July 15 th 2004). The contours are drawn differing in a factor of 10°' 25 ppmv. 
Figure 5. SABER vibrational temperatures of the main CO 2 t >2 vibrational states contributing 
to the 15-/xm channel typical of: a) mid-latitudes; b) polar summer; c) polar winter. The profiles 
are zonal means for 10° latitude boxes at the average latitude shown. The solid line is the kinetic 
temperature. 

Figure 6. Reported values for the relaxation rate of CC^t^) by molecular nitrogen. The 
solid line corresponds to the value adopted in vl.07 and the dashed lines indicate the uncertainty 
range assumed in this work. 

Figure 7. As in Fig. 6 but for the relaxation rate of C02(f2) by molecular oxygen. 
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Figure 8. Errors in SABER T*, caused by the uncertainty in the CC^t^) quenching rate by 
N 2 and O 2 for March 18 th (equinox, left) and July 15 th (solstice, right), 2004 (top: &co 2 - a ir x0.75; 
bottom: &cc>2-airXL25). These errors correspond to the temperature maps of Fig. 3. Contours 
are plotted every 0.5 K. 

Figure 9. Reported values for the relaxation rate of CO2(01 1 0) by atomic oxygen. 

Figure 10. As in Fig. 8 but for the uncertainty in CC^t^) quenching rate by atomic oxygen 
(top: fcco 2 -o x0.5; bottom: fcco2-oxl.5). Contours are plotted every IK. 

Figure 11. Effect on the C02(f2) vibrational temperatures of reducing the quenching rate of 
CO 2 with atomic oxygen by a factor of 2 (solid: nominal; dashed: perturbed) for typical mid- 
latitude (a) and polar summer (b) conditions. Differences in the left panels correspond to the 
perturbed minus the nominal vibrational temperature profiles for each level. Lines with symbols 
correspond to the vibrational levels indicated in the legend. The kinetic temperature is also 
shown in the left panels. 

Figure 12. Effect of the ^-quanta V-V exchange on the C02(f2) vibrational temperatures 
(solid: nominal k vv ; dotted: fc w =1.2xl0 -11 cm 3 s -1 ; dashed: k vv —2.4x 10 -11 cm 3 s _1 ) for typical 
polar summer conditions. Kinetic temperature is also shown. 

Figure 13. As in Fig. 8 but for the uncertainty in CO 2 V-V ^-quanta exchange rate (top: 
fcw=1.2xl0 -11 cm 3 s -1 ; bottom: fc vv =2.4x 10 -11 cm 3 s _1 ). Contours are plotted every IK. 
Figure 14. CO 2 vibrational temperatures for the main emitting levels at 15 /m n for a scan 
taken on January 15 th 2003. Labels as described in Fig. 5. 

Figure 15. Error in the kinetic temperature of the scan in Fig. 14 due to minimum (short 
dash) and maximum (long dash) values assumed for the quenching rate of 002 (^ 2 ) by N 2 and 
O 2 (left), by O (center) and by the V-V ^-quanta exchange rate (right). Note the different scale 
of the x-axes. 
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Figure 16 . Left: Average temperature for the ascending (solid) and descending (dash-dotted) 
SABER scans around the equator during the north viewing phase for July 15 th 2004. Center: 
Difference (ascending-descending) divided by 2 for the nominal case (solid), the minimum (short 
dashes) and the maximum (long dashes) values of the collisional rates. Right: Overall error in 
te difference divided by 2. 


Figure IT. SABER kinetic temperature (solid) and 002 (^ 2 ) vibrational temperatures (labels 
as described in Fig. 5) averaged for measurements taken around 80°N during January 15 th 2004. 
The anomalously high stratopause provokes the CO 2 levels to depart from LTE at unusual low 
altitudes. 


Table 1 . Collisional processes for the CO 2 vibrational levels relevant to this work, correspond- 


ing nominal collisional rates and associated uncertainty factors. 


No. Process® 

Colhsional rate b 

Uncertainty 

factor 

1: CO|(u 2 )+N 2 CO|(u2-1)+N 2 

^C02-airl— A^[7x10 17 VT+ 

+6.7xlO- 10 exp (-83.8/T 1 / 3 )] 

0.75-1.25 

2: C02(u2)+02 ^ C0 2 (i*2 — 1)+02 

^C02-air2— A^[7x10 17 VT+ 

+1 .0 x 10- 9 exp(-83.8 /T 1 / 3 )] 

0.75-1.25 

3: C0|(u2)+0( 3 P) # C0|(d 2 -1)+0( 3 P) 

^co2— o ~ 3 . 5 x 10 13 VT+ 

+2.3xl0- 9 exp(-76.75/T 1 / 3 ) 

0.5-1. 5 

4: CO*(v 2 )+COg(i/ 2 ) # COi(u 2 -l)+C(y 2 (^+l) 

fcw=1.2xl0 -11 if i—j &; v' 2 —0 

=2.4xl0 -11 if V 2 —V 2 +I—I 

See text 


a i and j denote different CO 2 isotopes. 


b Rates in the forward sense of the process in cm 3 s 1 . A V2= i — 1 and A V2 ^i— 3^, where g V2 is 
the degeneracy of level v^. 
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Table 2. Uncertainties in retrieved kinetic temperature (in K) due to the specified uncertainty 
in the collisional rates as in Table 1 for mid-latitude conditions. The columns labeled £ typ’ are 
the typical errors due to each rate, calculated with the quadratic average of the errors in both 
senses. TOTAL is the combined error due to the three rate constants in which: max - (estimated 
using Jsq of— air) ^col-o> an d below 88 km and k™ x above) and max + (estimated using 
^co 2 — ain &co 2 -o> an d Kv X below 88 km and k™ n above) are the extreme errors, and ‘typ’ is the 
root-square-sum of the typical errors due to each rate. 


z &C02— air &C02-0 & TOTAL 


[km] 

-25% 

25% 


typ 

-50% 

50% 


typ 

1.2 b 

2.4 b 


typ 

max 

max + 


typ 

104 

0.0 

0.0 

± 

0.0 

11 

-3.8 

± 

8.3 

0.1 

-0.2 

± 

0.1 

-3.9 

n 

± 

8.3 

100 

0.0 

0.0 

± 

o 

o 

4.8 

-1.7 

± 

3.6 

0.1 

-0.1 

± 

0.1 

-1.8 

4.7 

± 

3.6 

96 

0.0 

0.0 

± 

o 

o 

1.9 

-0.6 

± 

1.4 

0.1 

-0.1 

± 

0.1 

-0.6 

2.0 

± 

1.4 

92 

0.0 

0.0 

± 

o 

o 

1.7 

-0.5 

± 

1.3 

0.1 

-0.1 

± 

0.1 

-0.5 

1.8 

± 

1.3 

88 

0.0 

0.0 

± 

o 

o 

1.7 

-0.5 

± 

1.3 

0.0 

0.0 

± 

o 

o 

-0.5 

1.8 

± 

1.3 

84 

0.1 

-0.1 

± 

0.1 

1.0 

-0.5 

± 

O 

bo 

-0.3 

0.1 

± 

0.2 

-0.8 

1.3 

± 

0.8 

80 

0.1 

-0.1 

± 

0.1 

0.3 

-0.2 

± 

0.3 

-0.5 

0.1 

± 

0.3 

-0.8 

0.5 

± 

0.4 

76 

0.1 

0.0 

± 

0.1 

0.0 

0.0 

± 

O 

O 

-0.1 

0.0 

± 

0.1 

-0.3 

0.1 

± 

0.1 


a The Tfc bias produced by a 75% decrease in kco 2 -o, as suggested by laboratory measure- 


ments, is 0.7, 1.9, 6 and 16 K at 80, 84, 92, and 100km respectively (see text). 
b Rates in 10 -11 cm 3 s -1 . 
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Table 3. Uncertainties in retrieved kinetic temperature (in K) due to the specified uncertainty 


in the collisional rate for polar summer conditions. We used the minimum values of the rate 
constants below 88km and the maximum values above to estimate max - , and vice versa to 
estimate max + . 


Z &C02— air ^C02— O ^vv TOTAL 


[km] 

-25% 

25% 


typ 

-50% 

50% 

typ 

1.2 b 

2.4 b 


typ 

max 

max + 

typ 

104 

1.8 

-0.9 

± 

1.4 

43 

-16 

± 32 

4.6 

-2.2 

± 

3.6 

-18 

55 

± 33 

100 

0.9 

-0.5 

± 

0.7 

23 

-9.8 

± 18 

2.5 

-1.3 

± 

2.0 

-11 

25 

± 18 

96 

0.9 

-0.5 

± 

0.8 

14 

-6.3 

± 11 

2.7 

-1.3 

± 

2.1 

-7.6 

15 

± 11 

92 

0.8 

-0.5 

± 

0.7 

6.6 

-2.6 

± 5.0 

2.7 

-1.3 

± 

2.1 

-4.0 

8.1 

± 5.5 

88 

-0.4 

0.1 

± 

0.3 

1.0 

0.3 

± 0.7 

0.6 

0.2 

± 

0.5 

-0.1 

0.7 

± 0.9 

84 

-2.3 

1.6 

± 

2.0 

-2.5 

1.5 

± 2.1 

-7.1 

3.4 

± 

5.6 

-8.5 

6.0 

± 6.2 

80 

-1.1 

0.6 

± 

0.9 

-0.3 

-0.1 

± 0.2 

-4.1 

1.2 

± 

3.0 

-5.5 

1.8 

± 3.2 

76 

-0.2 

0.0 

± 

0.1 

-0.1 

-0.2 

± 0.1 

-0.7 

0.3 

± 

0.6 

-1.5 

0.4 

± 0.6 


a The Tfc bias produced by a 75% decrease in kco 2 -o, as suggested by laboratory measure- 


ments, is 0.1, -5, 14 and 51 K at 80, 84, 92, and 100km respectively (see text). 
b Rates in 10 -11 cm 3 s -1 . 


Table 4. Uncertainties in retrieved kinetic temperature (in K) due to the specified uncertainty 
in the collisional rate for polar winter conditions. The set of rate constant values producing the 
maximum errors max - and max + are as in mid-latitudes (see Table 2). 


Z &C02— air ^C02— O ^vv TOTAL 


[km] 

-25% 

25% 


typ 

-50% 

50% 


typ 

1.2 b 

2.4 b 


typ 

max 

max + 


typ 

104 

0.0 

0.0 

± 

0.0 

10 

-3.2 

± 

7.4 

-0.1 

-1.0 

± 

0.7 

-3.2 

n 

± 

7.4 

100 

-0.1 

0.0 

± 

0.1 

6.1 

-2.3 

± 

4.6 

-0.2 

-1.0 

± 

O 

bo 

-2.7 

6.7 

± 

4.7 

96 

-0.1 

0.0 

± 

o 

o 

2.8 

-0.9 

± 

2.1 

-0.8 

-1.1 

± 

0.9 

-1.9 

2.8 

± 

2.3 

92 

-0.1 

0.1 

± 

0.1 

0.5 

0.0 

± 

0.4 

-0.2 

-0.4 

± 

0.4 

-0.3 

0.5 

± 

0.5 

88 

0.0 

0.1 

± 

0.1 

-0.2 

0.3 

± 

0.2 

-0.4 

-0.2 

± 

0.3 

-0.7 

0.5 

± 

0.4 

84 

-0.1 

-0.1 

± 

0.1 

1.5 

-0.8 

± 

1.2 

-0.5 

-0.1 

± 

0.4 

-0.8 

1.8 

± 

1.2 

80 

0.5 

-0.4 

± 

0.4 

0.7 

-0.4 

± 

0.6 

0.1 

-0.2 

± 

0.2 

-0.9 

1.3 

± 

0.8 

76 

0.4 

-0.3 

± 

0.4 

0.1 

-0.1 

± 

0.1 

0.4 

-0.4 

± 

0.4 

-0.7 

0.9 

± 

0.5 


a The Tfc bias produced by a 75% decrease in kco 2 -o, as suggested by laboratory measure- 


ments, is 1.3, 3, 1.5 and 14 K at 80, 84, 92, and 100km respectively (see text). 
b Rates in 10 -11 cm 3 s -1 . 
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